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ABSTRACT 

The fraction of high-redshift sources which are multiply- imaged by intervening galaxies 
is strongly dependent on the cosmological constant, and so can be a useful probe of the 
cosmological model. However its power is limited by various systematic (and random) 
uncertainties in the calculation of lensing probabilities, one of the most important of 
which is the dynamical normalisation of elliptical galaxies. Assuming ellipticals' mass 
distributions can be modelled as isothermal spheres, the mass normalisation depends 
on: the velocity anisotropy; the luminosity density; the core radius; and the area over 
which the velocity dispersion is measured. The differences in the lensing probability 
and optical depth produced by using the correct normalisation can be comparable to 
the differences between even the most extreme cosmological models. The existing data 
is not sufhcient to determine the correct normalisation with enough certainty to allow 
lensing statistics to be used to their full potential. However, as the correct lensing 
probability is almost certainly higher than is usually assumed, upper bounds on the 
cosmological constant are not weakened by these possibilities. 

Key words: gravitational lensing ~ galaxies: statistics ~ galaxies: kinematics and 
dynamics - galaxies: structure - cosmology: miscellaneous. 



1 INTRODUCTION 

The fraction of high-redshift quasars which are multiply- 
imaged due to gravitational lensing is determined mainly 
by the cosmological model and the population of potential 
lenses. Many of the early investigations into the statistics of 
quasar lensing (e.g. Press & Gunn 1973; Turner, Ostriker & 
Gott 1984; Kochanek & Blandford 1987; Fukugita & Turner 
1991; Maoz & Rix 1993) focussed on the deflector popula- 
tion, but more recent studies have emphasised the cosmolog- 
ical possibilities. Specifically, Turner (1990) and Fukugita, 
Futumase & Kasai (1990) found that the lensing probabil- 
ity increases very rapidly with the (normalised) cosmological 
constant, O.Ao, but depends only weakly on the (normalised) 
matter density, Slmo • One of the most stringent upper limits 
that can be placed on the value of the cosmological con- 
stant is due to the low number of lenses detected - both 
Kochanek (1996b) and Falco, Kochanek & Munoz (1998) 
found that JIaq ^ 0.65, with 95 per cent confidence. These 
results are only marginally consistent with a number of in- 
dependent cosmological measurements, such as high-redshift 
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supernova observations (e.g. Schmidt et al. 1998; Perlmut- 
ter et al. 1999) and cosmic microwave background measure- 
ments (e.g. Lineweaver 1998; Efstathiou et al. 1999), which 
imply that = 0.7 ± 0.2, again at the 95 per cent con- 
fidence level. Further, the low density implied by cluster 
observations (e.g. Bahcall, Fan & Gen 1997), combined with 
the inflationary requirement of a flat universe (e.g. Guth 
1981; Kolb & Turner 1989) also imply a high value of ^Aq- 
It is thus very important to accurately assess both the ran- 
dom and systematic uncertainties on the lensing constraints. 

Some of the random uncertainties in the lens statistics 
are being steadily reduced as new surveys better constrain 
the deflector and source populations and more lenses are 
discovered. If ellipticals do completely dominate the lens- 
ing probability (e.g. Turner et al. 1984; Kochanek 1996b), 
any improvements in the knowledge of the number density 
of galaxies must be accompanied by accurate type informa- 
tion. Both the the Sloan Digital Sky Survey (e.g. Szalay 
1998; Loveday & Pier 1998) and the 2 degree Field galaxy 
redshift survey (e.g. CoUess 1999; Folkes et al. 1999) should 
decrease the errors on the type-specific luminosity functions 
by up to an order of magnitude. These two projects will 
also greatly reduce the uncertainties in the quasar luminos- 
ity function, as well as yielding more lensed quasars than 
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are known to date (e.g. Loveday & Pier 1998; Boyle et al. 
1999a,b; Mortlock & Webster 2000b). 

However, comparable progress in reducing the various 
systematic uncertainties is unlikely to be as rapid or as 
certain. Firstly, dust in the lensing galaxies can obscure 
multiply-imaged quasars from lens surveys. There is some 
evidence that it is an unimportant effect (e.g. Kochanek et 
al. 1999; Faico et al. 1999), but it is also possible that it 
dominates the statistics (Malhotra, Rhoads & Turner 1997) . 
There are a number of studies of dust in local galaxies, but it 
is reddening in high-redshift galaxies that is more important 
to lensing statistics. The only measurements of obscuration 
in such galaxies comes from lensed quasars, as the colours of 
the various images of the one source can be compared. Us- 
ing this technique Falco et al. (1999) found that ellipticals 
with redshifts of up to ~ 1 have minimal dust content (the 
difference in the extinctions between different lines-of-sight 
being only AE{B - V) ~ 0.2 mag). 

Another potential limitation on the accuracy of lens 
statistics is uncertainty in the mass evolution of galaxies. 
Keeton, Kochanek & Falco (1998) and Kochanek et al. 
(2000) have used lens galaxies to measure the fundamental 
plane (e.g. Dressier et al. 1997) of field ellipticals at moder- 
ate redshifts, but very little could be inferred about the mass 
evolution of the population. Assuming the present-day pop- 
ulation of ellipticals formed from the mergers of spirals (or 
other smaller halos), the high-redshift deflector population 
should consist of a larger number of less massive objects. If 
the total mass in halos is conserved the lensing optical depth 
is independent of the evolution, but the average image sep- 
aration is decreased (e.g. Rix et al. 1994). Mao & Kochanek 
(1994) used this fact to show that the known quasar lenses 
were best explained if there was little or no evolution in the 
elliptical population to redshifts of order unity. Thus a non- 
evolving population of elliptical galaxies is adopted here. 

The mass profile of the defiectors has a greater im- 
pact on the frequency of multiply-imaged sources, as well 
as the resultant image configurations. Constant mass-to- 
light ratio de Vaucouleurs (1948) models of ellipticals can be 
matched to either the galaxy dynamics (Kormendy & Djor- 
govski 1989; van der Marel 1991) or lens statistics (Maoz 
& Rix 1993; Kochanek 1996b), but not to both simulta- 
neously. The mass-to-light ratios required to reproduce the 
observed image separations are approximately double those 
suggested by dynamical arguments. This implies that el- 
lipticals are dominated by dark matter halos, which might 
be expected to follow the Navarro, Frenk & White (1996, 
1997) mass profile inferred from A^'-body simulations of cold 
dark matter-dominated galaxy formation. However, in the 
inner regions - which determine the strong lensing prop- 
erties - this profile is only marginally steeper than the de 
Vaucouleurs (1948) model, and is inconsistent with lensing 
observations for the same reasons. The inferred dark matter 
halos can be modelled as isothermal spheres (e.g. Binney & 
Tremaine 1987), which are consistent with both dynamical 
considerations (e.g. Kormendy & Richstone 1995) and lens- 
ing data (e.g. Kochanek 1993, 1996b). There is, however, 
some uncertainty as to the correct mass normalisation of 
this model (characterised by a velocity dispersion, (Joo), due 
to both its relationship with the observed line-of-sight ve- 
locity dispersion, (T||, and the possibility of a finite core. The 
lensing cross-section of an isothermal galaxy is proportional 



to (Turner et al. 1984), so even small variations in the 
normalisation are important. 

The surface brightness of ellipticals is fiatter than an 
isothermal profile near the centre (and steeper at large radii, 
although this is less relevant), which results in higher ob- 
served velocity dispersions for a given mass distribution than 
constant mass-to-light ratio models (e.g. Binney & Tremaine 
1987). Gott (1977) used the conversion = (3/2)^/V|| to 
account for the extended nature of the dark matter halo, 
but more recent calculations imply that the correct scaling 
is much closer to unity (e.g. Kochanek 1993, 1994). This re- 
sult is supported by lens statistics (e.g. Kochanek 1996b), 
but the uncertainties are quite large. 

If ellipticals do have finite cores (within which the den- 
sity is roughly constant), the the maximum deflection angle 
is reduced, making them less effective lenses (e.g. Hinshaw 
& Krauss 1987; Blandford & Kochanek 1987). However, the 
mass normalisation is increased for a given observed (T||, 
as the central potential well of the galaxies are shallower; 
this tends to increase their lensing effectiveness (Kochanek 
1996a,b). Further, non-singular lenses tend to produce more 
highly magnified images (e.g. Blandford & Kochanek 1987), 
resulting in an increased lensing probability due to magni- 
fication bias (Turner 1980). The qualitative arguments are 
quite clear, but the relative importance of the various effects, 
and their overall impact on lens statistics, are not. 

In Section ^ the normalisation and resultant scalings 
of the mass distribution are derived, and the effects these 
have on the optical depth and lensing probability of ellip- 
tical galaxies are discussed in Section ^ The conclusions 
reached on the effect of core radii on lens statistics are then 
summarised in Section ^ 



2 ELLIPTICAL GALAXIES 

A simple model for the population of elliptical galaxies (Sec- 
tion 2.1) is adopted here, in which individual galaxies are 
assumed to be spherically symmetric objects, completely de- 
fined by their radial mass distribution (Section 2.2), radial 



luminosity density (Section 2.3), and dynamics. Given that 



the line-of-sight velocity dispersion in the central regions is 
usually used for the dynamical normalisation of ellipticals, 
this observable must be related to the model parameters for 
self-consistency (Section 2.4). 



2.1 Population 

If elliptical galaxies follow a Schechter (1976) luminosity 
function, and obey the Faber-Jackson (1976) relationship, 
their local co-moving number density is given by 



7(l + a)-l 



exp 



(1) 



where o-|| is the observed line-of-sight velocity dispersion, 
a = -1.07 ± 0.05 and n, = (0.0019 ± 0.003)/i^ Mpc"^ (Ef- 
stathiou, Ellis & Peterson 1988). Here Ho = WOh km s"^ 
Mpc~^ is Hubble's constant, and cr, = 225 ± 20 km s~^ 
and 7 — 3.7 ± 1 (de Vaucouleurs & Olson 1982). There is 
the possibility of systematic uncertainties in equation (0) for 
low o-|| (e.g. Folkes et al. 1999), but the larger galaxies dom- 
inate the strong lensing by ellipticals. Further, any scatter 
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in the Faber- Jackson (1976) relation effectively increases ct, 
by an amount comparable to the scatter in (T|| (Kochanek 
1994). 

Under the assumption that the galaxy population is 
non-evolving (See Section |l|), the differential number of 
galaxies at redshift z and and with velocity dispersion (ji 



dz d(7|| 



dVo drig 
dz; d(T|| 



(2) 



where dVo/dz is the co- moving volume element at redshift 
z. Its full cosmological dependence is rather complex (e.g. 
Carroll, Press & Turner 1992; Kayser, Helbig & Schramm 
1997) and so only three simple, limiting cases are used here: 
f^mo = 1 and SIaq = (the Einstein-de Sitter model); f2mo = 
and JIaq = (the empty Milne model); and 57mo = and 
f^Ao = 1 (a cosmological constant-dominated flat model) . In 
these models the volume element becomes 



dVo 

d« 



(3) 



3 4(z-t-2^T+7+2) 
' (1+1)572 ' 



if r2mo = 1 and Qao = 0, 



47r 



47r I 



,Ho) 4{z + l)a 



if Vin 



and Qaq = 0, 



Ho) 



if f2mo = and SIaq = 1- 



2.2 Mass distribution 



Despite its unbounded total mass, the non-singular isother- 
mal sphere is consistent with the dynamics of elliptical 
galaxies (and their lensing properties). The Hinshaw & 
Krauss (1987) model has a mass density given by 



pM{r) 



(4) 



27rG (r2 + rl) ' 

where is the core radius and a^c is the line-of-sight ve- 
locity dispersion away from the core for a constant mass-to- 
light ratio galaxy. The integrated mass is given by 



M(< r) 



G 



arctan 



(a 



and the projected surface density by 

,2 -I 

Em(-R) 



2G Vi?2 + ri ' 
which can also be integrated analytically to give 



M(< R) = ^ (/r^ 



(5) 



(6) 



(7) 



As discussed in Section 2.3, the luminosity density of 



many ellipticals appears to be effectively singular, but such 
observations cannot directly constrain the mass distribution. 
Nonetheless, there are at least two strong arguments to sug- 
gest that Tc is small as well. Firstly, most lensed quasars have 
even numbers of images (e.g. Keeton & Kochanek 1996), 
which implies that the galaxies' mass distributions are very 
nearly singular (e.g. Wallington & Narayan 1993; Kassi- 
ola & Kovner 1993). Secondly, dynamical modelling, com- 
bined with high-resolution Hubble Space Telescope imaging 
of nearby ellipticals, reveals that a number have large black 
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Figure 1. The projected luminosity density of galaxies described 
by a de Vaucouleurs (1948) law (solid line) and a Hernquist (1990) 
profile (dashed line), scaled by their values at their effective or 
half-light radius, Rg. 



holes at their centres, and so are formally infinitely dense 
there (e.g. Kormendy et al. 1996, 1997). If the core radius 
is non-zero, its scaling with a is important (Kochanek 1991; 
Section and is taken to be 

r^-r.,(^y\ (8) 
where ttc = 4 ± 1 (e.g. Fukugita & Turner 1991). 

2.3 Light distribution 

Two models of the surface brightness of ellipticals are used 
here; both are compatible with observations, and the differ- 
ence between the results of the two models is an indication 
of the uncertainty of this calculation. 

The first is a de Vaucouleurs (1948) profile, given by 



EL(i?) = 296.7-^2 exp 



-7.67 ( — 



1/4 



(9) 



where _Rg is the effective or half-light radius of the galaxy 
and L is its luminosity. This is shown in Fig. |l| as the solid 
line. The luminosity density is given by an Abel integral 
(Binney & Tremaine 1987) as 



PL(r) = 



dSr 



dR ^ 



:dR. 



(10) 



It must be computed numerically for most r, but can be 
approximated by pL(r) ~ 1096.6 ///(tt^E^) {r/Rg)'^''^ for 
r-*0 (Young 1976). 

The second model is based on a Hernquist (1990) pro- 
file, which was developed as an analytical approximation 
to the de Vaucouleurs (1948) profile, but fits the data just 
as well in its own right. With Rn ~ 0.55-Rg,^ the surface 
brightness is 

t The definition Rjj — 0.45-Rg is sometimes used (e.g. Kochanek 
1996b), but it is more relevant for the dynamics of constant mass- 
to-light models than for luminosity profiles alone. 
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^ ru^ L [2 + {R/Rur][MR/Rn)-3] 



27ri?2 



[1 - {R/Ray 



where 



arccosh{l / x) 

1, 

arccos( 1 /x) 



if X <1, 



if X = 1, 



if X > 1. 



(12) 



This is shown as the dashed line in Fig. 0, and is qualitatively 
similar to the de Vaucouleurs (1948) profile for 0.1-Rg ^ R ^ 
10-Rg. The fraction of the flux in the discrepant regions is 
only a few per cent. The resultant luminosity density is given 
by 

= 27ri?3 r/Rn{l + r/Rn)^- ^^^^ 

The spatial scale of both distributions is determined 
by the effective radius, which, like the core radius of the 
mass distribution, is assumed to increase with the velocity 
dispersion of the galaxy as 

i?g = 7?g, (^)"', (14) 

where _Rg* = (4 ± l)/i kpc and = 4 ± 1 (Kormendy & 
Djorgovski 1989). 



2.4 Dynamical normalisation 

The pop ulation of galaxies is given in terms of (7| | in Sec- 



tion 2.1, but the mass distribution, and hence lensing prop- 
erties of individual galaxies are determined by a^o (Sec- 
tion gj). For the central dispersion decreases 
with increasing core radius, but also depends on both the 
stellar dynamics within the galaxy and the luminosity pro- 
file. Radial orbits result in higher dispersions, as the stars fall 
through the core of the galaxy, and more extended luminos- 
ity profiles also result in faster central stellar motions. The 
galax y m odel is s pecified by Pm(t') and Ph{r), given in Sec- 



tions 2.2 



and 



2.3, respectively, and the (assumed constant) 



velocity anisotropy, (3a- ■ This is defined as /9o- = 1 
where ae, (which are equal in a non-rotating system) 
and ar are the angular and radial components, respec- 
tively, of the velocity dispersion tensor of the luminous mat- 
ter. Both theoretical and observational results suggest that 
0.0 ^ /3a ^ 0.5 for ellipticals (Binney & Tremaine 1987; 
van der Marel 1991; Kochanek 1994), but a broader range 
of values is explored here. 

Under the above assumptions, the Jeans equation re- 
duces to (e.g. Binney & Tremaine 1987) 



ar (r) 



2 dcr^ ^ 1 



dpL ^ 2/3a 
(Tr(r) dr Phil") dr r 



GM{< r) 



which can be integrated to give 

-2/3^ 



(Jr{r) 



PL(r) 



GM(< r')pL(r')r 



/X /2{/3„-l) 



1 1/2 



dr 



(15) 



(16) 



This can be found analytically in the regions where pm is 
purely isothermal and Pl can be approximated by a global 
power law, given by p-L{r) oc [and hence El oc R~^^~^\ 



provided that ^ > 1]. The Hernquist (1990) model ap- 
proaches ^ = 1 for small r and ^ = 4 for large r; the de 
Vaucouleurs (1948) law has ^ = 3/4 for small radii. The 
radial dispersion in these regimes is then 



ar{r) 



1/2 



(17) 



In the outer regions the core radius has no effect, but for 
r ^ rc the lowered central density is more important than 
both PL and the velocity anisotropy (assuming f3a ^ 0.5). If 
Per — 1, the orbits are predominantly radial, and luminous 
material 'plunges' through the central regions of the galaxy, 
greatly increasing ar- 

The line-of-sight velocity dispersion at a given position 
on the galaxy is given by a projection integral (e.g. Binney 
& Tremaine 1987) as 



l(^) = 



EL(i?) 



{l-l3aRyr^)r 



i?2 



PL(r)ar(r) dr 



1/2 



(18) 



The power law approximations discussed above are also valid 
here; insertion of equation mTl) into the integral yields 



o-||(7?) = ao 



e - (C - l)/3. 



1/2 



(19) 



In the special case of ^ = 3, the standard result (Gott 1977), 
valid for all /3ct, that cr^^{R) = (2/3)^'''^ aoo, is recovered. It 
is also clear that the line-of-sight dispersion is lower in re- 
gions with a steeper luminosity density. Fig. ^ shows (T||/(Too 
as a function of R for the de Vaucouleurs (1948) and Hern- 
quist (1990) models. For both singular models and those 
with Per — 1, the central dispersion results is much higher 
than the dispersions at large radii. Conversely, the central 
dispersion decreases with core radius, as the galaxy's gravi- 
tational well is shallower, and bound orbits must be slower. 
Again the core radius is more important than the anisotropy 
of luminosity profile, provided only that Pa ^ 0.5. 

A real velocity dispersion is measured from a spectrum 
taken over a finite region of the galaxy, in the form of ei- 
ther a linear slit or, more commonly now, an optical fibre. 
Assuming a circular aperture of projected radius R{ centred 
on the galaxy, the observed line-of-sight velocity dispersion 
is 



o-||(< R{) = 



J^' 27ri?EL(7?)cr||(7£) dR 



/g^' 2tvRT,l{R) dR 



(20) 



Atmospheric seeing would have the effect of smoothing the 
function CT||(i?), but its importance is minimal as this tends 
to be a slowly varying function for R^ Rg. For a fixed angu- 
lar size, 6i, the aperture will vary from galaxy to galaxy as 
Rf = dA{0,z)6! ~ czOi/Ho, where z ^ 1 is the galaxy's 
redshift, and (Ia its the angular diameter distance. Only 
nearby galaxies are practical targets for dynamical stud- 
ies - for instance the samples of van der Marel (1991) and 
Lauer et al. (1995) contain galaxies with 0.002 ^ z 0.02 
and 0.0006 ^ z: ^ 0.06, respectively. Hence an aperture a 
few arcsec in diameter implies 0.01 kpc ^ Rt 1 kpc, as 
compared to core radii of ~ 0.1 kpc and effective radii of 
~ 4 kpc. Figs. I and § show the normalisation for singu- 
lar models as a function of Rf and Pa, respectively. The 
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Figure 2. The line-of-sight velocity dispersion, f7||, of the luminous matter in an elliptical galaxy with luminosity densities described by 
a de Vaucouleurs (1948) law in (a) and (b) and a Hernquist (1990) model in (c) and (d). In (a) and (c) the mass distribution is a singular 
isothermal sphere; in (b) and (d) it has a core radius of O.lRg, where Rg is the effective radius of the galaxy. The four lines in each panel 
represent different values of the velocity anisotropy: f3cr = 1-0 (solid lines); /3a- = 0.5 (dashed lines); /3a- = 0.0 (dot-dashed lines); and 
/3a = —0.5 (dotted lines). The regions where cr|| is flat are those where the luminosity density can be approximated as a global power 
law. 



normalisation increases with Rf as the observed dispersion 
is less dependent on the extreme orbital speeds near the 
core. Conversely, the normalisation decreases with (3a, as the 
maximum orbital speeds are higher for radial orbits. Fig. ^ 
shows croo/o"||(< -Rf) as a function of core radius. The gen- 
eral trend is that the observed velocity dispersion decreases 
(and so the dynamical normalisation increases) with core ra- 
dius and also with Rf. It is also clear that the quantitative 
behaviour varies comparably with the dynamical model and 
luminosity density assumed. For smaller apertures this de- 
pendence is more extreme, as only the core region is directly 
observed. 

Calculation of cr||(< Rt) is computationally expensive, 
especially for the de Vaucouleurs (1948) model; fortunately 
its dependence on core radius is very nearly linear. For the 
subsequent lens calculations, the normalisation used is 

aoo =<7|,(<i?t) (^A + B^^ , (21) 

where A and B independent of core radius. Kochanek 



(1996b) used A = 1 and 5 = 2 for galaxies with Rf — Rg, 
but this underestimates the normalisation by up to ~ 15 
per cent (See Fig. ^.), and hence underestimates the lensing 
probability by up to ~ 80 per cent. As defined above, A is 
determined purely by the singular models, for which the nor- 
malisation in plotted in Figs. ^ and ^ Also, A should match 
the power-law approximation given by equation (uM for very 
small apertures, and (Jcx)/o"||(< Rf) (3/2)^'^ ~ 1.225 
for J?f — > oo, as shown by Kochanek (1993). For interme- 
diate apertures, the variation between the results for the 
two luminosity profiles can be used as a guide to the uncer- 
tainty in the above results. Given this discrepancy, analyt- 
ical expressions for A and B were developed for both the 
de Vaucouleurs (1948) and Hernquist (1990) profiles, and 
these were used in the lensing calculations presented in Sec- 
tion ^. The analytical forms are given in Mortlock (1999), 
and, for all reasonable galactic models, agree with the nu- 
merical results to within a few per cent; the uncertainty in 
the normalisation caused by the ambiguity in the choice of 
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Figure 3. The dynamical normalisation, a^ , as a function of the size of the integration aperture, _Rf , for a singular galaxy with effective 
radius Rg and observed line-of-sight velocity dispersion (T||. The luminosity density is described by a de Vaucouleurs (1948) law in (a) 
and a Hernquist (1990) model in (b). Different values of the velocity anisotropy are denoted by the different line-styles: /3o- = 1-0 (solid 
lines); /Bcr = 0.5 (dashed lines); /Bcr = 0.0 (dot-dashed lines); and /3o- = —0.5 (dotted lines). 
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-0.5 
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Figure 4. The dynamical normalisation, aac, as a function of the velocity anisotropy, /So-, for a singular galaxy with effective radius 
Rg and observed line-of-sight velocity dispersion a^^. The luminosity density is described by a de Vaucouleurs (1948) law in (a) and a 
Hernquist (1990) model in (b). The lines represent different aperture sizes - from the bottom to the top, Rf/Rg takes on the values: 
0.0001, 0.001, 0.01, 0.1, 1.0 and 10.0. 



luminosity profile is considerably larger at around 10 per 
cent. 

Having found a relationsiiip between the isotliermal 
mass normalisation and observables for a given galaxy, tiiere 
still remains tiie question of how best to make the conversion 
from the observed Faber- Jackson (1976) relation. A 'natural' 
choice is to consider the aperture to be a specified fraction 
of the effective radius of each galaxy in the sample although 
this is unrealistic for large surveys undertaken using fibres of 
a fixed radius. (However, if the data is available, it is prefer- 
able to measure the velocity dispersion out to ~ i?g, as the 
results are more robust.) The other possibility is to treat 
measured dispersions as being averaged over a given physi- 
cal scale, independent of cr| | , but this breaks down for small 
galaxies as only the central regions are likely to be registered 
due to surface brightness considerations. One way to circum- 



vent these ambiguities is to fit the mass model directly from 
dynamical measurements on a galaxy- by- galaxy basis (e.g. 
van der Marel 1991; Kockanek 1994). This is not only more 
time-consuming, but also requires spatially resolved surface 
brightness and velocity dispersion measurements. 

The 'default' model (/3^ = 0; 7?f ~ Rg; = 0) has 
cToo/o"!! ~ 1.1, leading to a 50 per cent increase in the lens- 
ing probability, due to its cr^ dependence (e.g. Turner et al. 
1984; Kochanek 1994). For most of the other plausible sets 
of parameter values the differences are even greater. Simi- 
lar increases in the lensing probability can also result from 
the spread in the Faber- Jackson (1976) relation - a disper- 
sion of A (Too / (Too — 0.2 increases the lensing optical depth 
(Section |3.2| ) by a factor of ~ 1 -I- 6{Aaoo/o'ao)'^ — 1-25 
(Kochanek 1994). If a number of other common simplifi- 
cations used in lensing calculations are correct (e.g. that 
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Figure 5. The dynamical normalisation, (Too, as a function of core radius, r^, for a galaxy with effective radius Hg and observed line- 
of-sight velocity dispersion (T||. The luminosity density is described by a de Vaucouleurs (1948) law in (a) and a Hernquist (1990) model 
in (b). Different values of the velocity anisotropy are denoted by the different line-styles: = 1.0 (solid lines); P„ = 0.5 (dashed line); 
/3cr = 0.0 (dot-dashed lines); and /3<j = —0.5 (dotted lines). The two lines for each dynamical model are for different values of the aperture 
size: i?f = O.liJg; and R{ = Rg. The larger projected fibre radius results in less variation between dynamical models and higher values 
of (Too for Tc = 0, but lower values of CToo if the core radius is large. 



spiral galaxies are unimportant; that obscuration by dust 
is minimal; etc.), the dynamical normalisation is probably 
the greatest uncertainty in lens statistics at present. Fur- 
ther, whereas uncertainty in n* can be greatly decreased 
by simple counting procedures, accurate normalisation (or, 
equivalently, measurement of dng/d(Too) requires much more 
detailed data. 



3 LENSING PROBABILITY 

From the introduction of finite cores to lens galaxies by Hin- 
shaw & Krauss (1987) it was generally believed that the 
fraction of quasars which are multiply-imaged drops sharply 
with increasing core radius. However as shown by Kochanek 
(1996a,b) the requirements of self-consistent dynamical nor- 
malisation (Section 2.4) imply that the lensing probability 
declines only slowly, or possibly even increases with cor e ra - 
dius. Once the lens equation has been solved (Section 3.1), 
these possibilities are quantified in terms of both the optical 
depth (Section B.S) and a simple but self-consistent calcula- 



tion of lensing probability of a distant quasar (Section 3.3) 



3.1 The lens equation 

The lens equation relates the angular position (relative to 
the optical axis joining the observer and the centre of the 
lens) of the source, /3, to those of its image(s), 9. These 
quantities can only be related via dA{0,z^), dA(0, Zs) and 
dA{zd,Zs), the angular diameter distances from observer to 
deflector, observer to source, and deflector to source, respec- 
tively. In the flUed-beam approximation (Dyer & Roeder 
1972, 1973), the angular diameter distances in the cosmo- 
logical models described in Section 2.1 are given by 



dA{zi, Z2) 



(22) 



if Q,ma ~ and Q,Aa = 0, 



e Z2(z2 + 2)-zi{zi+2) 
Ho 2(zi+l)(z2+iy-' ' 



e Z2-ZI 
L Ho Z2 + 1 ■ 



if Qrng = and Qao ~ 1- 



For a thin lens with a cumulative surface mass distribu- 
tion M{< R), standard techniques (e.g. Schneider, Ehlers & 
Falco 1992) give the lens equation as 



(3^9- 



where 



M[< dA(0,^d) 



Tt9 d^(0, 2d)Scrit(.Zd,; 



Scrit {Zd: ^s) — 



dA(0, 



47rG dA(0, ^d) dA{zd, Zs 



(23) 



(24) 



is the surface density required for a rotationally-symmetric 
lens to be capable of forming multiple images (e.g. Subra- 
manian & Cowling 1986; Schneider et al. 1992). 

Using the mass distribution given in equation (Q), the 
lens equation becomes 



f3^9-9E 



where 9c = Vc/dAiQ, Zd) and 



= 47r 



if) 



dAjZd, Zs 
rfA(0, Zb] 



(25) 



(26) 



is the Einstein angle^ of the singular (rc = 0) lens. The dy- 
namical n orm alisation, (Too, is given in terms of observables 
in Section 2.4. 



t If the source, lens and observer are colinear, a circular image is 
be formed. The angular radius of this circle is the Einstein angle. 
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Equation ( |25| ) and 'a little algebra' (Hinshaw & Krauss 
1987) give 







2136' + [13' + 26e9c 



9l)e - 20e6Ic/3 



(27) 



where p = (36*1 + p' - 66Ie6Ic)/9 and q = /3(96l| + gfeeic - 
/32)/27. There are three solutions for S if < P"^, two if 
= and only one if > p^- Note that not all solutions 
need not correspond to image positions, as equation ( p7| ) is 
not strictly equivalent to equation (|2^). 

If 0c > 0e/2, equation ( |2^ ) can have only ever have 
one solution, as the central surface density of the lens is less 
than Ecrit [equation (^|)]. However, all but the most nearby 
galaxies have 6c < 0e/2 and so can form multiple images. 
In this case, expanding p and q, and disregarding multiple 
solutions that do not correspond to images, yields the result 
that a source is multiply- imaged is /3 < /3crit = /3-, where 



1/2 



If Tc 



(28) 

if 6*0 < 6^/2, 
if dc > Oe/2. 



0, the lens is always critical, and fUc: 



Be. 



cross-section for multiple-imaging is then simply 71(3^- 



The 
but 



the image positions and magnifications are needed for the 
more realistic lensing calculation (Section 3.S). If 6c <^ 6e 
and a"|| is fixed, /3crit decreases slowly with core radius, as 
there is little For larger core radii, however, 

(Too increases nearly linearly with rc, so, from equation (p6|), 
/3crit oc oc r^. Hence /3crit eventually increases with core 
radius. For large galaxies this only occurs for unrealistically 
large values of rc, but it can come into effect for very small 
galaxies if the fibre integration area is small (Section 3.2). 

If (0 <) /3 < /Scrit, the source is lensed, and the three 
images are located at 



HP) = -/3 + 2^cos 



+ 



2iTv 
~3~ 



(29) 



where i —1, 2 or 3, and two of the images are coincident if 
P = /Scrit- Note that this formula breaks down if the lens is 
singular, in which case only two images are formed. They 
are located at 6i — (3 + 9e and 62 — P — 9e- To within 10 
per cent, the angular separation between the outer pair of 
images is independent of source position, and is given by 
(Hinshaw & Krauss 1987) 



0, 



if Be > eE/2, 



AO : 



6Iev/1-26Ic/6Ie, if 6»c < fe/2. 



(30) 



This expression becomes exact in the singular case. 

If /3 > Peril, the source is unlensed, and there is only 
one image. For the range /3crit = P- < P ^ P+ the single 
image position is 



m ^lp+{q- V^^) +p(q- Vl^-P^) 



-1/3 



(31) 



If /3 > P+, equation (|2^) has three solutions, two of which 
do not represent image positions. The one solution that does 

is 



&iP) = -/3 + 2v^cos 



1 


f . ^ 




- arccos 




) 









(32) 



Again the singular case is much simpler: if (3 > 9e, the single 
image is located at 9 = f3 + Oe- 

The lens mapping (from 9 to f3, or image to source) 
changes the area subtended by an object; this manifests itself 
as an increase or decrease in flux, as surface brightness is 
conserved by gravitational lensing. The area change is given 
by the Jacobian of the mapping, and so the magnification is 
the reciprocal of this. For spherically-symmetric lenses (See 
Schneider et al. 1992.) the magnification of a point source is 
given by 

Pdp 



m(0) 

Equation (|25|) then gives 



(33) 



V0^ + Oi 



+ 9e9c - 

e^V9^TW 



(34) 



9e9? + (02 _ eE9c) ve' + ffi 

This expression breaks down for for small if ^c < 6'e/2 
but a Taylor expansion about 9 = gives 



A9l/9l 



(35) 



(20c/6'e - 1)2 + (26Ic/6(e - 1)02/02 . 

for ^ 0E. If fc = 0, equation ( ^^ reduces to pi(0) = 
ie/(i01-0E)l. 

The total magnification of a source is simply the sum 
of the magnification of its images, and is thus 



4/3)]. 



(36) 



The flux ratio is less well-defined, due to potential presence 
of three images and observational effects (e.g. Kassiola & 
Kovner 1993). For three image configurations it is taken to 
be the flux ratio of the two brightest images. 

3.2 Optical depth 

The lensing optical depth, r, as introduced by Turner et al. 
(1984), is the fraction of the source plane within which the 
lens equation has multiple solutions. It is a useful estimate of 
the lensing probability that is independent of observational 
restrictions and the source luminosity function. In the fiUed- 
beam approximation, the contribution to the optical depth 
by any one lens galaxy, Tg, is the fraction of the sky covered 
by its cross-section. Hence Tg = 7r/3crit/(47r), with /3crit given 
in equation (^8|) . The optical depth is given by integrating Tg 
over the population of deflectors, under the assumption that 
the individual cross-sections do not overlap. For a source at 
redshift 2:=, 



dz, dcrii ^ " ^ 



(37) 



If the lenses are singular and Oac ~ <^\\ is assumed, Tg = 0e/4 
and the integrals can be performed analytically (e.g. Turner 
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Figure 6. The 'standard' lensing optical depth, rgjg, as a func- 
tion of source redshift, Zb- This results from a population of sin- 
gular isothermal spheres with = <'"||i ^^^d the default parame- 
ters described in Section |2| Three cosmological models are shown: 
f2mo = 1 and Oaq = (solid line); Qmg = and Qaq = (dashed 
line); and flmo = and flAg = 1 (dot-dashed line). 



et al. 1984; Schneider et al. 1992; Krauss & White 1992; 
Kochanek 1993) to give 



TSIS 



(..) = 16.V(|^)^(^)V(l + a+i 



(38) 



X < 



, if f2mo = 1 and S7ao = 0, 

2(z^+4zf + 10z^-|-12za-|-6) ln(za + l)-3za(za-|-2)(z^+2za+2) 



8zJ(z3+2)2 



if Qr, 



■ and Qa 



^Zb, if Qmo = and ^Aq = 1, 

where equations (^, (^, (^), ( p^ ) and ( ^ ) have been used. 
The prefactor is the same as the _F-parameter introduced 
by Turner et aL (1984) to scale the lensing effectiveness 
of a galaxy populati on; f or the elliptical galaxy population 
described in Section 



2.1 



F ^ 0.008 ± 0.002. Fig. g shows 

TsisC^s) for the three cosmological models, showing the ex- 
pected increase with redshift as well as the strong depen- 
dence on JIaq- 

The optical depth (normalised to the above analytic 
cases) is shown as a function of core radius in Fig. |^. In- 
dependent of cosmology, the linear conversion aoa — 
results in t/tsis — 1.5; if the conversion is aoo — l-2cr||, 
then r/rsis ~ 2.1, and the optical depth is doubled. The 
lower set of lines in Fig. ^ show the marked decrease in r 
with core radius that was first demonstrated by Hinshaw 
& Krauss (1987). The upper set of lines use the normali- 
sation described in Section ^, and show not only the higher 
optical depth for singular models, but a very different depen- 
dence on core radius. For large rc, the effects of the normal- 
isation becomes very important, and 'small' galaxies (with 
cr|| ^ (T,) dominate the optical depth as their measured dis- 
persion is purely that in the under-dense core. The up-turn 
in r is not necessarily realistic, and occurs for core radii 
that are greater than observed - Kochanek (1996b) found 




C*/ g* 



Figure 7. The lensing optical depth, t, of a source at = 2, 
scaled by the 'standard' lensing optical depth for a population of 
singular isothermal spheres with ~ "^H- Three cosmological 
models are shown: = 1 and ^Aq = (solid lines); = 

and Haq = (dashed lines); and r2nio = a^nd Haq = 1 (dot- 



dashed lines). The lower set of lines are for (Too = 
of the other parameters are taken to be; /3ct = 0; 
and Rf = Rg. 



(T| I . The values 



: 4; Ug : 



Tc < 0.08i?g at 95 per cent confidence. For moderate core 
radii the normalisation is still important - in the fimo ~ 
and JIao = 1 cosmology, r/rsis ^ 1 for all values of rc. It is 
also potentially important that a finite core radius actually 
enhances the cosmological dependence of r, as most of the 
lenses are at high- redshift if flAg > (Kochanek 1992). 

The dependence of r on the various galaxy population 
parameters is shown in Mortlock (1999), but there is little 
difference between the resultant plots and those shown in 
Fig. ^ for the quasar lensing probability, and so the optical 
depth plots are omitted here. 



3.3 Quasar lensing probability 

The optical depth is a useful measure of lensing likelihood, 
but it cannot be directly compared to measured lensing fre- 
quencies. Hence the probability that a quasar is observed 
to be lensed, pq, must be calculated. In the case of the 
generic lens survey, all the quasars in a parent survey are 
re-examined for secondary lensed images. Under these con- 
ditions, pq is the fraction of all redshift Zq quasars of magni- 
tude mq (as measured in the parent survey) that would be 
revealed as lenses if examined with the resolution and sensi- 
tivity of the secondary search. The resolution limit implies 
a minimum image separation, A^min, and the depth of the 
follow-up observations leads to a maximum magnitude dif- 
ference between the primary and secondary images, Ammax. 
In general there is also a maximum image separation, A6'max, 
defined by the extent of the search for companion images. 

Most lensed sources are magnified by a factor of 2 or 
more, so the number of lenses at magnitude mq is deter- 
mined by the quasar number counts at least one magnitude 
fainter. The quasar luminosity function is so steep that this 
magnification bias (Turner 1980) can double or triple the es- 
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timated lensing probability. The quasar luminosity function 
is taken to satisfy 

j2 , 



Az dm 



(39) 



at all redshifts, where mqo = 19.0 ± 0.2 is the magnitude of 
the break in the number counts, aq = 0.9±0.1 is the bright- 
end slope, I3q = 0.3±0.1 is the faint-end slope (Boyle, Shanks 
& Peterson 1988; Kochanek 1996b), and the normalisation 
unimportant. 

The probability a quasar is lensed by a particular galaxy 
is given by integrating over the source position as (e.g. 
Kochanek 1995, 1996b) 



Xf"" 2^/3S(/3) 



Pq,! 



m=mq + 5/2 log[fitot (/3)] 



d/3 



4tv 



dzn drrin 



(40) 



where S{P) is the selection function. It can be approximated 

by 

S{I3) = H[Am„,ax - Am(/?)] 

X H[Ae{p) - Ae,„i„] H[A6in.ax - Ae{(3)] , (41) 

where H{x) is the Heavyside step function. Hence S{l3) = 1 
if the images of a source at position /3 satisfy the resolu- 
tion and sensitivity limits, and is zero otherwise. The use of 
A'tot(/3) [given in equation (^ij)] in the calculation of the mag- 
nification bias implies that the all the images of the source 
were unresolved in the parent survey]^, and the pre-factor 
converts the expression from the surface density of lenses to 
the fraction of quasars which are lensed. 

Integrating pq^g over the deflector population yields (c.f. 
Kochanek 1996b) 



Pq = 



dVo dng 

— — Pq,gda||d.g 



(42) 



This is somewhat simpler for the standard singular lens mod- 
els (e.g. Kochanek 1993), but no closed form expression for 
Pq is available unless a simpler form of the quasar luminosity 
function is used as well. 

Fig. ^ shows the probability that an niq — 19 quasar 
is observed to be lensed by an elliptical galaxy as a func- 
tion of source redshift. The survey parameters chosen are 

A^niin = 1 arcsec, A^niax ^ oJ^ and Ammax = 3. These val- 
ues affect the overall probability considerably (e.g. Schneider 
et al. 1992), but do not strongly influence the relative de- 
pendence on the other parameters. The lens probability is 
consistently a factor of ~ 2 higher than the optical depth 
shown in Fig. |^. Although some lenses are lost due to the an- 
gular separation cut-off, many more are magnified into the 
survey. As with the optical depth, pq is strongly-dependent 
on the cosmological model - the similarity of Figs, q and ^ 
demonstrates the generic nature of the cosmological depen- 
dence. 



§ This is almost always a valid assumption in the case of galactic 
lenses, but can result in a serious over-estimate of the magni- 
fication bias for more massive deflectors (Mortlock & Webster 
2000a). 

^ Galactic lenses are incapable of producing image separations 
of more than a few arcsec, so ASmax is usually unimportant. 




Figure 8. The probability that a rriq = 19 quasar at redshift Zq 
is observed to be multiply-imaged, assuming that all secondary 
images brighter than m = 22 and separated by > 1 arcsec from 
the primary are found. This results from a population of singular 
isothermal spheres with Coo = and the default population 
parameters described in Section ^. Three cosmological models are 
shown: Quig = 1 and = (solid line); flmg = and Caq = 
(dashed line); and Quig = and = 1 (dot-dashed line). 



Fig. ^ shows the dependence of Pq on various galaxy 
population parameters. Firstly, Fig. |9| (a) shows the varia- 
tion of lensing probability with core radius, in analogy with 
the optical depth dependence shown in Fig. |^. In the case of 
the unnormalised model (the single solid line), the high-rc* 
cut-off is even more pronounced than for r, despite the ef- 
fects of the magnification bias. This occurs as the image sep- 
aration drops below A^min quite quickly. For the correctly 
normalised models (the set of three lines), the enhancement 
in Pq is even greater than that for the optical depth (Fig. 
This comes about both due to the magnification bias and 
the increased deflection angles of the smaller deflectors. The 
increase in AO for a given is another important effect of 
the dynamical normalisation. 

Fig. ^ (b) shows that the velocity anisotropy, f3^ is not 
as important as the core radius in lens statistics. Whilst the 
difference between = —1 and f3cr ~ 1 can be a factor of 
several, most results sug gest that Per — for ellipticals, as 



discussed in Section 2.4. The lensing probability decreases 
with Pa as the observed dispersion within a fixed mass dis- 
tribution increases with the dominance of radial orbits. 

The variation of core radius and effective radius with 
velocity dispersion is usually unimportant in lensing calcu- 
lations, as illustrated by the flat parts of the curves in Fig. ^ 
(c) and (d) as well as Krauss & White (1992). However, 
in some situations, the values of Uc and Ug [as defined in 
equations (^) and (^^, respectively] can be important (e.g. 
Kochanek 1991). The sharp increases in pq seen in Fig. B (c) 
and (d) occur when \ug — Uc\ ^ 3, and the smaller galax- 
ies have core radii comparable to both their effective radius 
and to the scale over which the dispersion is measured. The 
massive increase in pq is probably unrealistic, but a weaker 
form of the effect will occur. 

The greatest assumption in these calculations is in- 
volved with the choice of R{, the scale over which the line- 
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Figure 9. The probability that a niq 



19, 



of a population of singular lenses with aoo 



■ 2 quasar is Icnscd with > 1 arcsec. It is scaled by the 'standard' lensing probability 
II . Three cosmological models are shown in each panel: fimo = 1 and J^Aq = (solid 



lines); Omo = and Oaq = (dashed lines); and Qmo = and Caq = 1 (dot-dashed lines). The lower line in (a) is for Ooa = ! the two 
sets of lines in panels (b), (c), (d), (e) and (f) represent rc* = (single line as there is no cosmological dependence) and rc» = O.liJg* 
(three distinct lines). The different panels show the dependence of pq on: scale core radius, rc*, in (a); velocity anisotropy, Pa, in (b); 
core radius scaling exponent, Uc, in (c); effective radius scaling, Mg, in (d); and the aperture size, ijf, in (e) and (f). In (e) ijf = ijg, 
scaling with galaxy size; in (f) R{ = i?g, and is constant. The defaults values used are: jSa = 0; Mc = 4; Ug = 4; and Rf = Rg. 
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of-sight velocity dispersion is measured. Fig. |^ (e) and (f) 
show Pq/pq,sis as a function of the typical scale of Rt; in 
the former panel it scales with the effective radius of the 
galaxies, whereas in the latter it does not vary with the 
properties of the galaxy in question. If ellipticals are singu- 
lar, a given observed dispersion results in a higher optical 
depth if i?f ^ 7?g, as the orbital speeds are much lower in 
the outer regions. Conversely, if ellipticals have significant 
core radii, the optical depth is highest for small apertures, 
in which case the measured dispersion is only a fraction of 
the dynamical normalisation. These arguments are true ir- 
respective of how Ri scales with _Rg, as can be seen from the 
similarity of panels (e) and (f) in Fig. |^. 

3.3.1 Cosmological implications 

The above formulation for pq was extended to arbitrary cos- 
mologies, the relevant distance and volume element formulae 
for which are given in e.g. Carroll et al. (1992). Fig. shows 
the interdependence of pq on the cosmological model and rc* , 
both with and without the self-consistent dynamical nor- 
malisation. In the models with Qaq = shown in (a), the 
core radius is considerably more important than the value 
of ilmo , as expected. Further, the difference between the 
normalised and unnormalised models is greater than that 
between any flAg = cosmological model. The spatially- 
flat models shown in (b) have a much stronger cosmological 
dependence, but even then the core radius becomes more 
important as it approaches the effective radius of the lens 
galaxies. For more realistic values of r^*, a slightly larger 
cosmological constant is permitted than if rc* — 0, but not 
relative to the a^o ~ o"|| models. For instance, if a given data- 
set implied an upper limit of Haq — 0.7 for unnormalised 
singular lenses, the limit becomes lower (S^Ao — 0.5) with 
aoo — 1-2(T|| singular galaxies. Even if rc* = 0.1-Rg, the up- 
per limit on the cosmological constant would still be ~ 0.6. 
The application of self-consistent dynamics shows that the 
standard models with aoo ~ cr|| and no core radius provide 
the weakest (i.e. very conservative) upper limits on flAg- 



4 CONCLUSIONS 

If the mass distribution of elliptical galaxies is essentially 
isothermal, one of the biggest uncertainties in their effect 
as gravitational lenses is their dynamical normalisation, Uoo- 
The mass scale depends on their internal dynamics (given by 
(T|| and /3o-), the luminosity profile and the aperture used to 
calibrate the Faber- Jackson (1976) relation. A range of non- 
rotating, spherical galaxy models reveals that aoo — 1-1(T|| 
for singular models, but that 1 ^ (Too/o"|| ^ 2 if ellipticals 
have significant core radii. 

The dynamical normalisation can have a strong effect 
on both the lensing optical depth and the more correct lens- 
ing probability. In the case of the singular lens model, both 
are increased (by up to a factor of 2) by the application of 
the correct normalisation, irrespective of the cosmological 
model. The effect of the normalisation is even greater with 
the presence of a core radius. Both the optical depth and 
lensing probability increase with very large core radii, purely 
due to these self-consistency requirements. For a given ob- 
served Faber- Jackson (1976) relation, the optical depth can 





= i-n. 



"'0 "0 

Figure 10. Contour plots showing the probability that a ruq = 
19, Zq = 2 quasar is lensed with A6 > 1 arcsec, as a function of 
Qmo ^rid the canonical core radius, rc*. In (a) there is no cos- 
mological constant (i.e. Qaq =0); in (b) the universe is spatially 
fiat (i.e. Qmg + ^Aq = l)- The contours are spaced logarithmi- 
cally — three per decade. In (a) those for the unnormalised lens 
model (dashed lines) range from pq = 2 X 10~^ to pq = 5 X 10~^; 
those for the correctly-normalised model (solid lines) range from 



Pq 



5 X 10 * to Pq = 5 X 10 In (b) those for the unnor- 



malised lens model (dashed lines) range from Pq = 2 X 10"^ to 
Pq = 2 X 10"'^; those for the correctly-normalised model (solid 
lines) range from pq = 5 X 10~^ to pq = 5 X lO"'^. 



vary by a factor of several with both the size of the disper- 
sion aperture and the velocity anisotropy. 

If lensing statistics are treated primarily as a cosmo- 
logical probe, the above uncertainties place limits on the 
accuracy of any cosmological inferences. However, because 
the correct normalisation almost certainly increases the cal- 
culated lensing probability in a given cosmological model, it 
strengthens arguments against a high cosmological constant. 
The weakest limits on Qaq obtained with the correct normal- 
isation and an arbitrarily large core radius are comparable 
to the limits obtained for the standard singular models with 
aoo — The dynamical conversions must be determined 
more precisely, theoretically and especially observationally, 
to make full use of gravitational lensing statistics. 
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